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Is there a link between the neutron-star spin and the 
frequency of the kilohertz quasi-periodic oscillations? 
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ABSTRACT 

There is a general consensus that the frequencies of the kilohertz Quasi-Periodic Oscil- 
lations (kHz QPOs) in neutron-star low-mass X-ray binaries are directly linked to the 
spin of the neutron star. The root of this idea is the apparent clustering of the ratio of 
the frequency difference of the kHz QPOs and the neutron-star spin frequency, Av/v s , 
at around 0.5 and 1 in ten systems for which these two quantities have been measured. 
Here we reexamine all available data of sources for which there exist measurements 
of two simultaneous kHz QPOs and spin frequencies, and we advance the possibility 
that Av and v s are not related to each other. We discuss ways in which this possibility 
could be tested with current and future observations. 
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1 INTRODUCTION 

Observations of neutron star low-mass X-ray binaries (NS 
LMXBs) with th e Rossi X-ray Timing Explorer (RXTE; 
iBradt et all 1 19931 ) have led to two important discoveries: 
Strong variability on millisecond timescales in the X-ray 
light curves of these systems, th e so-called kilohertz quas i- 
periodic oscillations (kHz QPOs: Ivan der Klis et al"1ll996al ). 
and pulsa tions during X-ray bursts , also known as burst os- 
cillations (jStrohmaver et al1ll996ah . 

The kHz QPOs are relatively narrow peaks in the power 
density spectrum of NS LMXBs that often appear in pairs, 
at frequencies v\ and V2 > v\ that change with time. These 
QPOs are thought to reflect motion of matter at the inner 
edge of an accretion disk around the neutron star. 

Burst oscillations are short-lived (r < 10s) , almost co- 
herent pulsations seen at the rise and tail of X-ray bursts 
in NS LMXBs. The frequency of these oscillations, Vb, in- 
creases in the tail of the bursts to an asymptotic value that 
is consistent with bein g the same in bursts se parated by 
more than a year time (|Strohmaver et al] [l998). This, and 
the fact that in the accretion-powered millisecond X-ray pul- 
sar (AMP) SAX J1808.4-3658 burst oscillations appear at 
the same frequency as the pulsations seen d uring persistent 
(non-burst) intervals l|Wiinands et al.ll2003l ). indicates that 
the frequency of these burst oscillations is equal to the spin 
frequency of the NS, v a . 

It is commonly accepted that the spin of the neutron 
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star is directly involved in the mechanism that produces 
the kHz QPOs. This consensus stems from the first detec- 
tion of kHz QPOs and burst oscillations in the same source, 
the LMXB 4U 1728-34, very early on in the RXTE mis- 
sion. While in different observations the kHz QPOs appeared 
at different frequencies, v\ in the range ~ 600 — 800 Hz, 
and V2 in the range ~ 500 — 1100 Hz, the frequency differ- 
ence of the QPOs, when both were present simultaneously, 
was consistent with being constant, Av — V2 — v\ » 363 
Hz, and also consistent with the oscilla tions seen during 
bursts in this source at v\, — 363 Hz (Strohmavcr et al 



1996b). This fitted with the suggestion IjStrohmaver et al 



1996c) that a beat mechanism with the neutron star spin 



was responsi ble for the kHz Q POs. Further results on other 
sources (e.g.. lFord et~ailll997l ) appeared to confirm this pic- 
ture. A detailed model, t he so nic-point model, proposed by 
iMiller. Lamb, fc Psaltis! (|l998l ) explained the observed rela- 
tion between the kHz QPOs and the neutron star spin in 
terms of a beat between material orbiting at the inner edge 
of the disk with the Keplerian frequency at that radius, and 
the spin of the NS. 

As soon as kHz QPOs were discovered in 4U 1636- 
53 l|van der Klis et al] [l996b) with a frequency difference 
of Av = 272 ±11 Hz, and burst oscillations at a frequency 
Vb = 581 Hz (|Zhang et al.|[l99^ ). it became apparent that 
in this source Av was inconsistent with being equal to Vb, 
but it was close to f&/2. This would have been the end of 
the sonic-point model, unless in 4U 1636-53 the 581 Hz fre- 
quency seen during X-ray bursts was the second harmonic 
of the NS spin frequency, Vb = 2 x i/ S) with v s — 290.5 
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Hz, e.g. if the pulsed radiation came from two antipodal 
poles on the NS. Although searches for a signal at half the 
burst oscillations frequency, the putative spin frequency of 
the neutron star ((Miller 1999), in the power spe ctrum of the 
bursts in 4U 1636-53 yielded no posit ive result l|Strohmaver] 
l200ll ; lstrohmaver fc Mark wardt 2002j), this option remained 
viable. 

When kHz QPOs and burst oscillations were discovered 
in more sources, it became apparent that there was a sys- 
tematic trend in how Af and i/j were related: For sources 
for which v\> < 400 Hz, Av ~ Vb, whereas for sources for 
which Vb jSj 400 Hz, Av ~ Vb/2. These two groups of sources 
were then called "slow" and "fast" rotators, respectively 
(|Muno et al.ll200ll ). 

Related to this, it is interesting to note that when plot- 
ted against each other, the frequencies of the kHz QPOs in 
19 diff erent sources all follow approximately the same re- 
lation jBelloni. Mendez. fc Homanll2005l . 120071 ; IZhang et all 
2006). This is a priori unexpected if in each source the fre- 
quency of the upper and lower kHz QPOs were related to 
the spin frequency as V2 = v\ + v s , given that v\ and V2 
span more or less the same frequency range in all sources of 
kHz QPOs, whereas the neutron stars in these systems have 
spin s that span a large range of frequencies, v s ~ 200 — 620 
Hz (|Chakrabartv et al J 12003 ). 

More kHz QPO data, and more precise QPO fre- 
quency measurements, showed that at least in some 
sources Av was not constant , but decreased a s the 
QPO frequencies increased (jvan der Klis et all 1 19971 ; 
iMendez et all Il999j^ and it was always s ignificantly lower 
than either ~b JlVlendez fc van der Klisl 19991 ) or ut /2 
II Mendez. van der Klis. fc van Paradiisl 119981). Mo difica- 
tions of the sonic-point model ( Lamb fc Mil ler 2001) could 
account for this difference, considering that the frequencies 
of the QPOs drift slightly when the material that produces 
the QPOs crosses the sonic point and falls onto the neutron 
star surface. 

Three other results raised more serious issues against 
the sonic-point beat-freq uency model, and eventually 
render ed it untenable: (i) Ijonker. Mendez fc van der Klisl 
(2002) found that in 4U 1636-53, when the frequency of the 
kHz QPOs decreases sufficiently, Av is significantly higher 
than Ub/2, which was difficult (if not impossible) to explain 
by the so nic-point model , even after the modifications intro- 
duced bv lLamb fc Miller! (l200ll); in the AMP SAX J 1808.4- 
3658, fii) Chakrabar tv et al.l (|2003l ) found that the frequency 
of burst oscillations is equ al to the NS spin frequency (see 
iMarkwardt fc Swan k 2003, for a similar result in another 
AMP, XTE J1814-338), while (hi) IWiinands et~ai1 (|2003l ) 
detected two simultaneous kHz QPOs with a frequency sep- 
aration Av ~ v a /2. If this is extended to other sources 
in which Av ~ Vb/2, it would also be true that for those 
vt = v s , and hence Av ~ v 3 /2. The sonic-point model could 
not explain this. 

But soon after the SAX J1808. 4-3658 results were pub- 
lished, a new model that reestablished a relation between 
the spin frequency of the NS an d the kHz QPO, the sonic- 
point and spin-res onance model |Lamb fc M iller 2003), was 
proposed (see also iLee et aL I l2004l ). In this model, there is 
a resonance in the accretion disk at the radial distance at 
which the Keplerian orbital frequency is equal to the neutron 
star spin frequency minus the vertical epicyclic frequency. 



This resonance could lead to either Av — v s or Av — v s /2 
depending on whether the disk flow at the resonance radius 
is smooth or clumped. In fact, the same source may in prin- 
ciple show both cases, but this has so far not been observed. 

Almost in parallel with some of these explanations, and 
as a result of some of the difficulties for beat-frequency mod- 
els mentioned above, a different class of models was pro- 
posed, in which the frequencies of the kHz QPOs were as- 
sociated to two of the three epicyclic freque ncies of general 
relati vity, or a combination of those (e.g., IStella fc Vietril 
1999). In these models, the NS spin frequency plays no role 
in setting up the frequencies of the kHz QPOs, except for the 
small corrections it introduces to the epicyclic frequencies. 
While these models reproduce qualitatively the trends seen 
in the da ta, and predicted other trends that were la ter on 
observed (|Migliari et al.ll2003l ; iBoutloukos et al.ll2006l '). they 
have problems to fit the data in detail. The main criticism to 
these models, however, has always been that they do not ex- 
lain the fac t that in several sources Av ~ v s or Av ~ v s /2 
Lambll2003l ). In other words, the criticism is that in these 
models the NS spin plays no role in the mechanism that 
produces the QPOs. 

Recently, lYin et all |2007l) compared the average fre- 
quency separation of the kHz QPOs, (Av), with v s in 
six systems in which these two quantities have been mea- 
sured. They suggest that, despite the low number of sources 
available for their analysis, (A^) depends weakly on v s , 
(Av) ~ -0.20v s + 390 Hz. 

In summary, the history of models of the kHz QPOs is 
a cycle of attempts to explain the phenomenon in relation 
to the spin of the NS; each time that a new observation 
raised an issue against one such model, a modification of 
that model, or a new model, was proposed that tried to 
reestablish the role of the NS spin in the mechanism that 
produces the kHz QPOs. 

After more than ten years from the discovery of the 
kHz QPOs, a critical assessment of the current paradigms 
is necessary. Here we review all the values of Av and v a 
available in the literature in order to compare the slow/fast 
rotator paradigm with other possibilities. We suggest that 
the data may in fact show that there is no relation between 
NS spin and kHz QPOs. Actually, the data appear to be 
consistent with a situation in which the average separation 
in frequency of the kHz QPOs (A^) is more or less constant, 
independent of the spin of the neutron star in the system. 
The division between "slow" and "fast" rotators may be an 
effect of the low number of sources for which two simulta- 
neous kHz QPOs and burst oscillations and/or pulsations 
in the persistent emission have been observed, and the fact 
that (Av) is independent of v s . 



2 DATA 

We use data from the literature. For the rest of the paper we 
assume that the frequency Vb of burst oscillations is equal to 
the spin frequency of the neutron star, v s (see JTJ). There are 
ten sources for which both Av and v s have been measured. 
Two of these sources are the AMPs SAX J 1808.4-3658 and 
XTE J1807-294, six of them are the atoll sources 4U 1608- 
52, 4U 1636-53, 4U 1702-43, 4U 1728-34, KS 173 1-260, 
and 4U 1915-05 (see lHasinger fc van der Klisl Il989l . for a 
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Table 1. Measurements of the frequency separation of the kHz 
QPOs, Av, and spin frequencies, v s , for the 2 AMP and 6 + 2 
atoll sources studied in this paper. 



Source 


Av range (Hz) v 


s (Hz) 


References 


Accretion-powered millisecond X-ray pulsars 


XTE J1807-294 


180-250 


191 


1 


SAX J1808.4-3658 


195 


101 


2 


Atoll sources 


4U 1608-52 


225-325 


619 


3 


4U 1636-53 


215-330 


581 


4 


4U 1702-43 


330 


330 


5 


4U 1728-34 


270-360 


363 


6 


4U 1731-260 


265 


524 


7 


ICR J17191-2821 


330 


294 


8 


SAX J1750.8-2900 


317 


600 


9 


4U 1915-05 


290-355 


270 


10 



REFE RENCES 
||2003|) 



(11 ILinares et al] 1120051): (21 IWiinands et al.l 



(31 IMendez et al.l ||19 99); Munoetal. 



(2001) 



(4) 



iDi Salvo et al.l [120031): Ijonker et al] ll2002h: IStrohmaver etaT 
||1998|1 ; (5) iMarkwardt et al] 1119991) : ( 61 IMendez fc van der Kli; 
||1999|1 : IStrohmaver et al.l dl996b() : (71 IWiinands fc van der Kli; 
|U997p ; ISmith, Morgan, fc Bradtl [1199711; (8) IMarkwardt et al 
1200711: iKlein-Wolt et al.l [1200711: (91 iKaaret et al.l l|2002l) ; (101 
iBoirin et al.l | |2000]1 ; iGallowav et al.l l l200ll) 



definition of the atoll class), and the other two are IGR 
J17191-2821 and SAX J1750. 8-2900 which most likely are 
also atoll sources. 

For each of these sources we give in Table [1] the spin 
frequency and the range of measurements of Av. Although 
there are too few sources to draw firm conclusions, from this 
Table it is apparent that the average of the Av range in the 
AMPs SAX J1808.4-3568 and XTE J1807-214 is somewhat 
lower than for the other eight sources. 

In Figure [T] we show in black the distribution of mea- 
surements of Av for the ten sources with known spin fre- 
quencies in Table [TJ and for the other sources of kHz QPOs 
for which the spin frequency is not known. The data for 
this Figure were taken from the papers in which the kHz 
QPOs were measured (see van der Klis 2006 for a complete 
reference list). We note that to measure both kHz QPOs 
significantly and calculate Av, in most cases power spec- 
tra had to be selected on the basis of some property of the 
source (e.g., intensity, colours, characteristic frequency of a 
low- or high-frequ ency timing feature, etc.; see for instance 
Ijonker et al. 2000) and averaged. In this process, part of the 
information of the real distribution of Av is lost. Therefore, 
the plots in Fig.[T]do not show the real Av distribution, but 
just the range of values observed and a rough idea of how of- 
ten a certain value of Av has been observed. This plot again 
suggests that the average Av of the AMP, (Av) ~ 210 Hz, 
is somewhat lower than the average Av of the other sources, 
(Av) ~ 300 Hz. 

Figure [2] shows the plot of Av/v s vs. v s for the ten 
sources m Tabled] For each source we plot all the individual 
Av measurements, ta ken from the refere nces listed in Table 
[1] divided by v s (see Ivan der Klisll2005l . for a similar plot). 
The open circles correspond to the two AMPs, and the filled 
circles are the other eight sources in Table [1] The dashed line 
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Figure 1. The distribution of measurements of Av in different 
groups of sources. As described in the text, because of the way 
the measurements were done, the histograms do not necessarily 
represent the true distribution of Av (see text for details). Upper 
panel: Atoll sources for which the spin frequency of the neutron 
star is known from burst oscillations. Second panel: Atoll sources 
for which the spin period of the neutron star is not known. Third 
panel: Z sources; spin frequencies are not known in this case. 
Lower panel: Accreting millisecond X-ray pulsars. The black his- 
togram shows the actual Av measurements; the gray histogram is 
the distribution of measurements mult i plied by th e factors close 
to 1.5 taken from Ivan Straaten et al.l l|2005l l and ILinares et al.l 
ll2005ll . 
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Figure 2. The ratio of individual measurements of Av divided 
by v s as a function of v s for the sources listed in Table [T] The 
open circles correspond to the two AMPs, and the filled circles 
are the other eight sources in Table [T] The dashed line is the step 
function, S{v a ) = 1 for v a < 400 Hz, S(v a ) = 0.5 for v a > 400. 



in this Figure is a step function, S(v a ) = 1 for v a ^ 400 Hz, 
S(v a ) = 0.5 for v s > 400. This Figure shows the fact that 
for "slow rotators", neutron stars with v a <400 Hz, the ratio 
Avjv a ~ 1, whereas for "fast rotators", neutron stars with 
v s >400 Hz, the ratio Av/v a ~ 0.5. Note that, as mentioned 
in SJT] for some sources Av is significantly different from v a 
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or v 3 /2, respectively, hence some of the individual ratios are 
signi ficantly different from 1 or 0.5, respectively. 



11999) 



Psaltis. Belloni. &; van der Klisl 



van Straaten et all 



van Straaten. van der Klis. fc Mended 



Reig, van Straaten & van der Klis 



2002), 



2003), 



2004), 



van Straaten, Wijnands & van der Klisl ([20051 ) . and 



Altamirano et al. 1 2007T ). have shown that there is a 



correlation between the frequency of the kHz QPOs and 
that of other low-frequency QPOs. More specifically, all 
these authors have shown that when plotted vs. the fre- 
quency of the upper kHz QPO, the frequency of the lower 
kHz QPO as well as the frequency of all low-frequency 
QPOs follow individual correlations that are consistent 
with being the same in five at oll sources, one Z source (see 
lHasinger fc van der Klisl Il989l . for the definition of the Z 
class), three low- lumino sity bursters, and two AMP s (see 
Ivan Straaten et al .1120051 , and I Altamirano et al.ll2007l . for an 
overview of these correlations). 

The AMPs SAX J1808.4-3568 and XTE J1807-214 
show relations between the frequencies of the low-frequency 
QPOs and V2 that are similar to those of the low-luminosity 
bursters and atoll sources. But the relations of SAX J1808.4- 
3568 and XTE J180 7-214 are shifted with respect to those 
of th e other sources (|van Straaten et al.ll2005l ; iLinares et al.l 
2005). The shiflQ is between the frequencies of the low- 
frequency QPOs and V2, and is best described as a mul- 
tiplication of U2 by a factor close to 1.5. The exact multi- 
plicative factors are 1.45 for SAX J1808.4-3568 and 1.59 for 
XTE J1807-214, respectively. While this factor applied to 
works for the low-frequency QPOs vs. V2 correlations, it 
does not work for the correlation betw een vi and V2 - Inter - 
estingly. Ivan Straaten et al.l (|2005h and lLinares et all (|2005l ) 
noticed that the v\ vs. V2 correlation in the AMPs and in 
the other sources could be reconciled if they also multiplied 
v\ by the same factor that they used to describe the shift of 
the V2 vs. low-frequency QPO correlations. (Notice that v\ 
was not used to derive that factor.) 

The nature of this shift is unclear, but taken at face 
value, a multiplicative factor applied both to v\ and V2 im- 
plies that the frequency difference Av = v?. — v\ must also be 
multiplied by this factor. The gray histogram in the lower 
panel of Figure [T] shows the distribution of Av (apart from 
the caveat described above in this section) for the two AMPs 
XTE J1807-294 and SAX J1808 4-3658 multiplied by the 
fact ors (close to 1 . 5) ta ken from Ivan Straaten et al.l (|2005l ) 
and lLinares et all (|2005l ). This multiplicative factors appear 
to bring the values of Av in these two AMP into the range of 
the values measured in all the other sources. The average Av 
of the combined sample (Atoll sources with spin frequency 
and the two AMPs multiplied by the factors close to 1.5) is 
308 Hz, and the standard deviation is 38 Hz. 

We note that the above procedure could imply a circular 
argument: Matching the v\ vs. V2 correlation of the AMPs 
and the other sources via a multiplicative factor means that 
also Av of the AMPs and of the other so urces would match. 
However, both Ivan Straaten et al.l (|2005l ) and lLinares et ail 



1 T his is usually called a shift since Ivan Straaten etafl (2005) 
and ILinares et al. measure a frequency shift in a log-log 

plot. 




Figure 3. The ratio of individual measurements of Av divided 
by v s as a function of v s for the sources listed in Table [T] but 
with Av of the AMP mu ltipli ed by the factors clo se to 1.5 from 
van Straaten et al ] l|2005l l and lLinares et al] {2005). The symbols 
and the dashed line are same as in Figure [2] The solid line cor- 
responds to a constant Av = 308 Hz. The dotted line shows the 
relation Av = -0.20v a + 390 Hz from lYin et al] [|2007r i. 



(2005) determined the shift factor on V2 using only the cor- 
relations between the low- frequency QPOs and V2, indepen- 
dently of v\. They then noted that they could also match 
the V\ vs. V2 correlations of the AMPs and the other sources 
if they applied the same factor (within errors) also to v\ 
(thence, there would be no free dom in choosing the shift fac- 
tor on v\\ see the description in ILinares et al Mil). Withal, 
the argument would indeed be circular if the shifts on v\ and 
on V2 turned out not to be the same. This conundrum may 
eventually be resolved when more shifts in other sources are 
observed. In the meantime, we caution the reader about the 
possible caveats in our procedur e of multiplying Av in the 
AM Ps by the fac t ors fo und by Ivan Straaten et all (|2005l ) 
and lLinares et aL I |2005h . 

In Figure [3] we show the plot of Av/v s for the ten 
sources in Tablefl] but now we have mult iplied the Av values 
of the AM Ps by the factors ta ken from Ivan Straaten et"all 
l|2005h and lLinares et al] (|2005l 'l. As in Figure H the dashed 
line is the step function S(v) (see above); the solid line cor- 
responds to a constant Av — 308 Hz. The dotted lin e show s 
the relation Av = -0.20v s + 390 Hz from lYin etHI l|2007h . 



3 DISCUSSION 

There is a general (but not universal) tendency to try and 
include the spin of the neutron star, f s ,asa key ingredient in 
models that explain the kHz QPOs. In this type of models, 
the spin is related to the difference between the frequencies 
of the kHz QPOs, Av — 1/2 — V\. This tendency persisted, 
despite the fact that several results seemed to contradict 
the predictions of these models. Amendments to the original 
ideas meant that models had to appeal to rather contrived 
geometries to explain the data, and rather artificial classes 
of sources had to be introduced to explain the diversity of 
the results (e.g., the division of sources with a different re- 
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Figure 4. Upper plot: The distribution of measurements of Av, 
of all types of sources in Fig[T] combined. The solid line repre- 
sents the best-fit Gaussian to this distribution, with a centroid of 
303.2±2.9 Hz and a standard deviation of 36.0±2.1 Hz (x 2 = 14.3 
for 12 degrees of freedom; 1-a errors shown). Lower plot: The dis- 
tribution of u s in twenty three sources, nine AMPs with pulsations 
in the per sistent emission and fourteen sources with burst oscilla- 
tions (see lYin et alj|2007l . and Table [T] for the list of sources and 
spin frequencies). 



lation between Av and v s into "slow" and "fast" rotators). 
Here we propose that the data are in fact consistent with a 
simpler picture in which the frequency separation between 
kHz QPOs, Av, is independent of v s , with (Av) more or less 
constant across sources. 

In Figure|4]we compare the distribution of the Av mea- 
surements in all types of sources for which two simultane- 
ous kHz QPOs have been detected (see Fig. [TJ, with the 
distribution of spin frequencies of twenty three sources for 
which pulsations in the persiste nt emission or burst oscil- 
lations have been measured f see lYin et al. I l2007l . and Table 
[1] for the list of sources and spin frequencies). From this 
Figure it is apparent that the distribution of Av measure- 
ments is much more concentrated than the distribution of 
spin frequencies. The distribution of Av measurements can 
be well described (\ 2 = 14.3 for 12 degrees of freedom) by 
a Gaussian with a mean value (Av) = 303.2 ± 2.9 Hz and 
a standard deviation cjav = 36.0 ± 2.1 Hz (1-a errors). A 
Kolmogorov-Smirnov (K-S) test yields a very low probabil- 
ity, P ~ 7 x 10~ 7 , that the two samples are drawn from the 
same parent population. According to the paradigm of slow 
and fast rotators (see SJTJ) , when v a >400 Hz (the exact value 
is not specified) Av should be compared to v a /2 instead of 
v s . After dividing by two the spin frequencies higher than 
400 Hz, a K-S test indicates that the distributions of v s and 
Av are marginally consistent with each other, with a K-S 
probability P ~ 1 x 10 -2 that the two are drawn from the 
same parent population. 

Using Av and v s for six of the a toll sources also in- 
cluded in our sample, lYin et ail l|2007h proposed that Av 
may be (weakly) related to v a in a way that is different 
than predicted by beat-frequency models. They found that 



(Av) ~ -0.20j/ s + 390 H20. In their analysis, lYin et all 
l|2007l ) treat the two AMPs differently, and do not include 
them in th eir fit. Following the results of lvan Straaten et al.l 
|2005l ) and lLinares et al l (|2005l ). here we do include the two 
AMPs in the analysis. We take a step further than lYin et al.l 
(2007). and we advance the hypothesis that in fact Av and 
v a are independent quantities. 

Our interpretation th at Av is independe n t of v a re- 
lies on the di scover y by Ivan Straaten et"al! (|2005h and 
iLinares et ail i|2005l) . who found that to reconcile the 
frequency-frequency correlations of the two AMPs SAX 
J1808.4-3568 and XTE J1807-214 with similar correla- 
tions in other sources, the frequencies of the kHz QPOs 
in the AMPs have to be multiplied by a factor of ~ 1.5 
(see fl2] for a discussion of p ossible caveats of th is). Both 
Ivan Straaten et al.l (120051 ) and lLinares et al.l l|2005l ) find that 
they can also reconcile the frequency-frequency correlations 
if they apply different multiplicative factors to the fre- 
qu ency of all variability com ponents, except V2- (As noted 
by Ivan Straaten et"aTI 120051 , a single multiplicative factor 
close to 1.5 applied both to v\ and V2 without changing the 
frequency of the other variability components is the sim- 
plest option.). In particular, they find that if V2 remains 
unchanged, vi of the AMPs SAX J1808.4-3568 and XTE 
J1807-214 has to be multiplied by a factor ~ 0.8 - 0.9 to 
match the v\ — V2 correlation defined by the atoll sources and 
low-luminosity bursters. The picture presented here does not 
change significantly if we only apply the ~ 0.8 — 0.9 factors 
to v\ and calculate new Av values. 

Although it is not the purpose of this paper to ex- 
plain the nature of these factors, here we provide some ideas 
about their possible origin. The usual suspect is the mag- 
netic field. A stronger field could prevent the inner edge 
of the disk from moving inward and, if the kHz QPOs are 
produced at that radius, a larger inner disk radius could 
imply lower kHz QPO frequencies, v\ and V2- If, on the 
contrary, the low-frequency variability is produced at larger 
radii, they would be less affected by the neutron-star mag- 
netic field, which could explain why a shift of the fre- 
quency of the low-frequency components is not required 
to match the frequency-frequency correlations. The prob- 
lem with this explanation is that other AMPs supposedly 
having relatively high magnetic fields, at least comparable 
to those in SAX J1808. 4-3568 and XTE J1807- 214, show 
no or very small shifts l|van Straaten et al.ll2005l ). Also, at 
least one other non-pulsating source shows shifts in the 
correlation, although sm aller than the ones in the AMPs: 
lAltamirano et all (|2005h find a shift of ~ 1.15 for the LMXB 



2 Notice that to get this result, lYin et at | [120071) used the in- 
dividual measurements of Av for each source to calculate (A^), 
and the spread of the Av values as a measure of the error. In the 
case of KS 1731-260 for which there is only one measurement of 
Av, they instead used the error of the measurement of Av. In 
the case of 4U 1702—43 all Av measurements were done over a 
narrow range of QPO frequencies, which may artificially reduce 
the spread. 

3 Kendall's r test (e.e. lPress et al.lll992l) appli ed to the d a ta pre - 
sented here, averaged in the same way as in lYin et al.l (120071) ■ 
yields a very low probability (between 1.7 and 1.9 a, depend- 
ing on whether we include the two AMPs are not) that there is 
correlation between (Av) and v s 
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4U 1820-30. The other sources showing significant shifts in 
the frequency-frequency correla tions are the AMPs XTE 
J0929-314, with u 3 = 185 Hz (iGallowav et~ai"1 120021 ) and 
a shift of 1.48 ± 0.11 (Ivan Straaten et al.ll2005f). and XTE 
J1814-338 with u a = 31 4 Hz jMarkwardt fc Swankll2003h . 
and a shift of 1.21 ± 0.09 |van Straaten et al.ll2005l . the shift 
here is marginally significant). 

A neutron-star mass differ ence could also e xplain these 
shifts. E.g., in the model of IStella fc Vietril |l999i), and 
in other models that explain the frequencies of the kHz 
QPOs in terms of epicyclic frequencies in general relativ- 
ity, the relation between u\ and vi depends explicitly on the 
neutron-star mass (see , e.g.. lStella. Vietri. fc Morsinkll 19991 ; 
iBoutloukos et al.ll2006l '): As noted bv lBelloni et all (^OOTt ). a 
multiplicative factor applied to the neutron-star mass trans- 
lates into the same mul tiplicative factor applied to both 
kHz QPOs (see eq. (4) in lStelfa et ailll999l '). If this (model- 
dependent) interpretation is correct, a factor ~ 1.5 in u\ and 
vi for the AMPs implies that the neutron stars in those sys- 
tems are ~ 1.5 times more massive than in the other atoll 
and Z sources and the low-luminosity bursters. It is generally 
assumed that, due to accretion, the neutron stars in LMXBs 
have masses larger than the canonical I.4M0 neutron star. 
If the ~ 1.5 factor is related to a difference in neutron-star 
mass, this would imply uncomfortably large masses for the 
neutron stars in the two AMPs. It is somewhat more dif- 
ficult to assess the effect of the mass of the neutron star 
on the low-frequency components, because it is not clear 
what frequency in the model represents the frequency of 
those components. If one of these low-frequency QPOs were 
due to Lense-Thirring precession ilStella fc Vietrilll998l). its 



frequency would be (see eq. (1) in IStella fc Vietril 1 19981 ) 
vlt oc IM~ x v\v s , where / and M are the moment of in- 
ertia and the mass of the neutron star, respectively, and as 
usual V2 is the frequency of the upper kHz QPO and u s is 
the spin frequency of the neutron star. If this identification 
is correct, the spin frequencies in table [1] a shift factor 1.5 in 
the mass and in the frequency of the upper kHz QPO, but 
no shift of the low-frequency QPOs imply that the moment 
of inertia must be a factor ~ 10 different among some of 
these sources. 

We know that Au is not the same in all sources 
jZhang et a l. 2006) and not even for the same source 
when more than one s ignificant detection is available 
(Ivan der Klis et al. 1 1 19971 ; iMendez et aT1ll999h . However, it 
is remarkable that, other than in the two AMPs, in all 
sources in which two simultaneous kHz QPOs have been 
detected, {Au) is approximately the sa me. After applying 
the m ultipl icative factors descr ibed in Ivan Straaten et"al] 
( 2005]) and iLinares et alj (|2005h , the same is true for the 
two AMPs (see for possible caveats) . This despite the 
fact that the measured spin frequencies span a fa ctor of 
more than 3. In the model of IStella fc Vietril (|l999f ). Au is 
equal to the radial epicyclic frequency, u r which, for the case 
of negligible eccentricity and a non-rotating neutron star is 
v r = (1 ~ 6GM/rc 2 ) 1/2 u^, with V4 = 1/ \2%){GM /r 3 ) 1/2 
the azimuthal frequency, identified with U2 in their model. 
(For neutron stars with spins smaller than ~ 600 Hz and 
masses in the range 1.4 — 2Mq, taking a typical range of U2 
frequencies, the radial epicyclic frequency is within » 15% 
of the value given by this formula.) One would then expect 
that on average Au would be the same for all neutron stars 



if they all have more or less the same mass and their upper 
kHz QPO spans more or less the same frequency range. 

If it is generally true that (Au) is more or less the 
same in all sources of kHz QPOs, the idea that there is 
no link between Au and u s could be tested in the case of the 
LMXB EXO 0748-676, which h as a spin frequency of 45 Hz 
|Villarreal fc Strohmarerll2004h . From the results in Figure 
[3] for EXO 0748-676 one expects Au « 300 Hz, whereas 
models that include a direct link between neutron-star spin 
frequency and frequencies of the kHz QPOs predict that for 
this source Au should be either 22.5 Hz or 45 Hz. In fact, in 
the context of "slow" and "fast" rotators, for EXO 0748-676 
Au is expected to be 45 Hz. (We cannot discard that if two 
simultaneous kHz QPOs are detected in EXO 0748-676 with 
Au ps 300 Hz, there would be attempts to modify existing 
models, or propose completely new ones, to explain Au/u s 
ratios that are an integer larger than 1.) Unfortunately, so 
far a single kHz QPO has b een observed from this source 
jHoman fc van der Klis]|200Ch . 

An equally interesting test of this idea would be to find 
a source with a spin frequency in the range u s ~ 350 — 500 
Hz for which no (or only a small) shi ft is required to fit 
the fr equency-frequency correlations of Ivan Straaten et all 
(2005). A case of more or less constant {Au) across differ- 
ent sources implies that Au/u s would be between 0.6 and 
0.8. Actually, there is a source that may be used for this 
in the near future: The AMP XTE J1751-305 h as a spin 
frequency u s = 435 Hz (|Markwardt et all l2002h , whereas 
Ivan Straaten et alj (|2005f ) find that a shift of only 1.12±0.03 
applied to U2 is required to match the frequency-frequency 
correlations. Unfortunately, so far there has been no detec- 
tion of two simultaneous kHz QPOs that would allow us to 
calculate Au in this source. 

Despite the fact that the data seem to suggest that 
{Au) is more or less the same in all sources of kHz QPOs 
(in the cases of the AMP SAX J1808.4-3658 and XTE 
J 1807-294 after applying a multiplicative factor deduced 
fro m the low-frequency QPO vs. U2 correlations; see above 
and Ivan Straaten et al. l l2005l andlL inares et al. 2005 for de- 
tails and possible caveats), we have no strong reason to dis- 
card the possibility that there are sources for which this is 
not the case (even after applying factors similar to those 
deduced in SAX J1808.4-3658 and XTE J1807-294). Our 
conjecture that there is no link between Au and u s would 
therefore not be weakened if a source with two simultaneous 
kHz QPOs is ever discovered, for which {Au) is not close to 
~ 300 Hz, as long as in such a source {Au) (after accounting 
for any possible shift factor as the ones in SAX J1808. 4-3658 
and X TE J1807-294; lvan Straaten et alj|2005l ; lLinares et all 
l2005h is different from u a and u s /2. 

From FigureEJit is apparent that in SAX J1808.4-3658 
and X TE J1807-294 the shifts o n the frequency of th e kHz 
QPOs |van Straaten et alj 120051 ; ILinares et al.ll2005h . and 
the idea that Au is either equal to u s or u s /2 are inconsistent 
with each other. It seems unlikely that this issue can be fully 
resolved as long as the nature of frequency shifts remains 
unexplained. While here we propose that the shifts imply 
that Au is not equal to u s or u s /2, we cannot completely 
discard that the shifts have a different explanation, and that 
in the two AMPs Au/u 3 is indeed close to either 1 or 0.5. 
We note, however, that even without taking the shifts into 
account, there is solid evidence that in several sources Au is 
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significantly different from i/ s or v a /2 (see JTJ). The question 
is how strong the evidence must be before the idea that Az/ 
and v a are directly linked is abandoned. 

To conclude, here we put forward the idea that the 
frequency difference of the kHz QPOs, Ais, in neutron- 
star low-mass X-ray binaries may not be related at all to 
the spin frequency, v 3 , of the neutron star. Beat-frequency 
mechanisms have been proposed not just in the context of 
the kHz QPOs; they we r e originally advan ced in the 1980s 
(|Alpar fc Shahaml 1 19851 ; lLamb et all Il985l ) to explain the 
low-frequency QPOs in these systems. We cannot rule out 
completely the hypothesis of a similar type of link between 
the kHz QPOs and the neutron-star spin, but this idea can 
in principle be tested and, if proven wrong, discarded. 
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